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Abstract 

We study statistical properties of the energy spectra of two-dimensional quasiperiodic tight-binding models. The 
multifractal nature of the eigenstates of these models is corroborated by the scaling of the participation numbers with the 
systems size. Hence one might have expected "critical" or "intermediate" statistics for the level-spacing distributions as 
observed at the metal-insulator transition in the three-dimensional Anderson model of disorder. However, our numerical 
results are in perfect agreement with the universal level-spacing distributions of the Gaussian orthogonal random matrix 
ensemble, including the distribution of spacings between second, third, and forth neighbour energy levels. 
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1. Introduction 

The anomalous transport properties of quasicrys- 
tals evince the peculiar electronic properties of these 
intermetallic alloys, which are reflected in the na- 
ture and the density of electronic states in quasiperi- 
odic systems. Tight-binding models on quasiperi- 
odic tilings serve as simple toy models of elec- 
trons moving in a quasiperiodic environment. While 
much is known rigorously for one-dimensional ape- 
riodic Schrodinger operators , quasiperiodic tight- 
binding models in the physically relevant dimensions 
two and three have mainly been investigated numer- 
ically, compare |^ and references therein for the few 
exceptions, which, however, yield only rather limited 
information about the overall structure of the spec- 
trum or the wave functions. 

At least in two dimensions, numerical results 
provide ample evidence that typical eigenstates of 
quasiperiodic tight-binding Hamiltonians are, as in 
the one-dimensional case, neither extended nor ex- 
ponentially localized; instead, the probability am- 
plitudes show a multifractal distribution Q. How- 
ever, the integrated density of states does not show 
the hierarchical plateau structure observed for one- 
dimensional substitution chains Q ; for pure hopping 
Hamiltonians it appears to be rather smooth, apart 



from a few gaps in the spectrum. This makes it pos- 
sible to employ spectral statistics Q as further tools 
to characterize the energy spectrum of quasiperiodic 
tight-binding models. 

2. Energy-level statistics 

The statistical analysis of the spectra of quan- 
tum systems has a rich history, and is at the cen- 
ter stage of current interest ||]. For disordered sys- 
tems, it is by now well known that in the metallic 
regime the level-spacing distribution ^0(5) of neigh- 
bouring eigenenergies is well described by P^'-'^(s) 
of the Gaussian orthogonal ensemble (GOE) of ran- 
dom matrix theory Q. Here, s denotes the en- 
ergy spacing in units of the mean level spacing, and 
Pq{s) gives the distribution of normalized gaps in 
the spectrum in the limit of infinite system size. In 
the insulating regime, the level-spacing distribution 
follows the Poisson law Pq{s) — exp(— s). At the 
metal-insulator transition, the level-spacing distribu- 
tion has been shown to follow yet another behavior 
which is attributed to the existence of a "critical en- 
semble" 

Recently, we have shown that the level statistics of 
tight-binding models defined on planar quasiperiodic 
tilings, notably the Ammann-Beenker tiling shown 
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Fig. 1. Top: Octagonal cluster of the Ammann-Beenker 
tiling with 833 vertices and exact Dg-symmetry around 
the central vertex. Bottom: A square-shaped cut with 
496 vertices without any exact spatial symmetries. 



in Fig. ^ is also given by the GOE predictions 
I IC I , contrary to previous statements in the hterature 
in, |2[ O. T he deviations from GOE behaviour ob- 
served in ]11|, could be understood by noting that 
the standard periodic approximants of the octagonal 
tiling are singular patches whose fourfold rotational 
symmetry is broken only weakly, see Fig. |^. This "al- 
most symmetry" of the patch does not lead to a block 
structure of the Hamiltonian (which would mean in- 
dependent spectra for each block), but nevertheless 
affects the level-spacing distribution. As shown in 
fldji , the level-spacing distribution Po{s) of generic 
patches of the tiling, or of the irreducible sectors of 
patches with exact symmetries, follows the GOE dis- 
tribution when applying a suitable unfolding proce- 
dure to correct for the fluctuations in the density 
of states. This result was further corroborated by 
higher correlations such as the Mehta-Dyson statis- 
tics A3 [|lOj and the number variance E2 Fur- 
thermore, we have shown that the level-spacing dis- 
tribution remains GOE if one considers small parts of 
the spectrum without any unfolding, and we demon- 
strated that the unfolding procedure has to be per- 
formed somewhat differently from the procedure in 



Fig. 2. Periodic approximant of the Ammann-Beenker 
tiling with 1393 vertices, displayed together with a copy 
of itself that has been rotated through 90 degrees. The 
"worms" where mismatches between the two rotated 
copies occur are shaded. 



disordered systems p5[ . 

In this paper, we numerically calculate not only 
the first, but also the second, the third, and the 
fourth neighbour level-spacing distributions P„(s), 
n = 0, ...,3, respectively. We show that the 
higher level-spacing distributions are also well de- 
scribed by the predictions of random matrix the- 
ory. The recently proposed "semi-Poisson ensem- 
ble" supposed to be valid for the multifrac- 
tal eigenstates of the metal-insulator transition in 
the three-dimensional Anderson model jl^, does 
not describe the higher level-spacing distributions 
of the quasiperiodic models. The same applies 
to other intermediate statistics obtained recently 
for an ensemble constructed by an analogy with a 
one-dimensional gas model of particles with nearest 
neighbour interactions V{x) = — log|x| ||l8|. This 
is intermediate between the Poisson and the GOE 
statistics, which can also be interpreted in terms 
of Boltzmann distributions of one-dimensional in- 
teracting gases with pair potentials V{x) = and 
V{x) = — log |a;|, but not restricted to nearest neigh- 
bours, respectively. Finally, we calculate the partic- 
ipation numbers in these tilings and show that their 
scaling behavior indicates multifractal behavior. 

3. Higher level-spacing distributions 

In addition to the nearest-neighbour level-spacing 
distribution Po{s) considered previously jlO[ [l5|] , 
we have now computed the distributions Pi{s), 
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Fig. 3. Level-spacing distribution up to the forth neigh- 
bour spacing for the largest sector of the octagonal patch 
with N = 157369 vertices. The thin lines represent the 
GOE predictions 

^2(5), and ^3(5) for the simple Hamiltonian 

N 

H^J2^^,\^){J\ (1) 

defined on the Ammann-Beenker or octagonal tiling 
as shown in Fig. |^. Here, each vertex i = 1, 2, . . . , TV 
of the tiling carries a state \i), and the hopping ele- 
ments tij = tji are chosen to be unity between neigh- 
bouring vertices, i.e. vertices connected by a bond, 
and zero otherwise. The spectral statistics are ob- 
tained by numerically diagonalizing the Hamiltonian 
(^) with free boundary conditions for large patches 
of the tiling, taking into account their symmetries. 
While the GOE predictions for these distributions 
are known in principle ||^, it is non-trivial to ob- 
tain numerical data for these functions. As we could 
not find sufhciently accurate results in the litera- 
ture, we calculated P^'^^{s), P^°^{s), P^°^[s), 
and P|^°^(s) using the approach of jl^ by nu- 
merically solving the appropriate differential equa- 
tions. This was done using Mathematica by per- 
forming a small-s expansion, followed by a numer- 
ical integration to larger values of s. The small-s 
behaviours of the GOE distributions are P^'~'^{s) — 
TT^s/e + 0(s3), pGOE(g) ^ ^4^4/270 + 0(56), 
PfOE(g-) ^ 7r8sVl764000 + O(si0), and P^°^{s) = 
TT^'^s^^ / 168781725000 + 0{s^^). In our calculations 
we evaluated the expansions beyond the 30th order. 

In Fig. ^, we show the results for the data of the 
largest patch considered in |l^, an octagonal patch 
as in Fig. ^but with 157369 vertices, employing the 
unfolding procedure described in |l^. Apparently, 
our data agree very well with the GOE distributions. 
We like to emphasize that the higher level-spacing 
distributions considered here, and in particular their 
small-s behaviour, are very sensitive to detect devia- 



Fig. 4. Comparison of the level-spacing distributions 
of Fig. ^ with the results for the semi-Poisson (dotted) 
and the one-dimensional gas model (dashed) ensembles. 
Note that Po{s) is identical for these two intermediate 
ensembles. 

tions from the GOE behaviour. However, we do not 
find deviations beyond the fluctuations expected for 
finite systems. Thus, the level-spacing distribution 
is the same as that of the Anderson model of local- 
ization in the metallic regime j^, ||. This is surpris- 
ing, because it is well-known that eigenstates in the 
metallic regime are extended, whereas there is ample 
evidence that the eigenstates on planar quasiperiodic 
tilings have a multifractal character. Thus one might 
speculate that the level statistics should have been 
similar to the intermediate statistics found jl^ at 
the metal-insulator transition in three dimensions, 
where eigenstates are multifractal as well. In Fig. ^, 
we compare our results to the level-spacing distri- 
butions of the intermediate semi-Poisson ensemble 
[ p^ [1^ and the one-dimensional gas ensemble with 
logarithmic interactions ||l^. These latter statistics 
clearly do not agree with our data. 

4. Scaling of the participation numbers 

We now turn our attention to the calculation of 
the participation numbers P{N^ E) in order to in- 
vestigate the properties of the states in the tiling di- 
rectly. For this purpose, we consider square-shaped 
patches as in Fig. [| with up to = 7785 vertices 
[ p^ and compute 

N 

p{N,Er^^Y.\'i^m\ (2) 

i=l 

where ^{E) — {'tpi{E),... ,iPn{E)) denotes a nor- 
malized eigenstate at energy E. Thus a state i^iE) 
that is completely localized at site J, i.e. ipiiE) = 
has P{N,E) — 1, whereas a fully extended state 
'iP^(E) = 1/y/N gives P{N,E) = A^. Because for 
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Fig. 5. Energy-averaged participation numbers versus 
system size for square-shaped patches (cp. Fig. |l|). The 
straight line is a least-squares fit of log P{N) = k log N+c 
to the data. 

the models under consideration the level-spacing dis- 
tribution does not appear to depend on the energy 
jlSt , we concentrate on the energy-averaged partic- 
ipation number P{N). For multifractal states, we 
expect to find a power-law behaviour P{N) ~ iV" 
as a function of systems size , between the result 
for localized (k — 0) and extended {k = 1) states. 
The result is shown in Fig. ^, which yields a value of 
K « 0.87 ± 0.05. This value, in agreement with pre- 
vious results for quasiperiodic tight-binding models 
1^, is clearly below the values k = 1 for extended 
states, albeit rather close, corroborating the mul- 
tifractal character of the eigenstates in our model, 
though with rather large characteristic exponents . 

5. Conclusions 

We considered the spectral properties of a tight- 
binding Hamiltonian defined on the Ammann- 
Beenker tiling. Although this system is not ener- 
getically disordered, the level-spacing distributions 
exhibit level repulsion and are in perfect agreement 
with the GOE predictions. This is similar as for the 
energetically disordered Anderson model of localiza- 
tion in the metallic regime and also for tight- 
binding models defined on certain random graphs 
||2l|] . Nevertheless, the scaling of the energy-averaged 
participation numbers confirms that the correspond- 
ing eigenstates are not extended. Their multifractal 
character is further corroborated by the system-size 
dependence of the participation numbers. At first 
glance, it appears surprising that multifractal states 
may give rise to GOE level statistics. However, it 
is known from studies of the metal-insulator transi- 
tion in more than two dimensions that the critical 
level statistics becomes similar to the Poissonian be- 
haviour for larger dimensions p2[ . Thus, it appears 
possible that the critical level statistics for the two- 
dimensional case considered here just coincides with 



the Gaussian orthogonal random matrix ensemble. 
This scenario is in perfect agreement with the results 
presented here. 
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